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As a direct source of information on chiral symmetry breaking within QCD, the sigma 
commutator is of considerable importance. Since hadron structure is a non-perturbative 
problem, numerical calculations on a space-time lattice are currently the only rigorous 
approach. With recent advances in the calculation of hadron masses within full QCD, 
it is of interest to see whether the sigma commutator can be calculated directly from 
the dependence of the nucleon mass on the input quark mass. We show that, provided 
the correct chiral behaviour of QCD is respected in the extrapolation to realistic quark 
masses, one can indeed obtain a fairly reliable determination of the sigma commutator 
using present lattice data. For two-flavour dynamical fermion QCD the sigma commutator 
lies between 45 and 55 MeV based on recent data from CP-PACS and UKQCD. 

1. WHAT IS THE SIGMA COMMUTATOR ? 

In the quest to understand hadron structure within QCD, small violations of funda- 
mental symmetries play a vital role. The sigma commutator, a^'- 

VN = l (N\ [Q a , [Qis, H]] | A) = m (N\ uu + dd \N) = m^- (1) 

(with Qi5 the two-flavour (i=l, 2, 3) axial charge) is an extremely important example of 
such a symmetry. 



2. PREVIOUS ATTEMPTS 



On cannot be accessed directly by experimental measurements. However, one can infer 
from world data a value of 45 ± 8 MeV This result has been under some scrutiny 
recently due to the progress in new determinations of the pion-nucleon scattering lengths 
PH and new phase shift analyses The full lattice QCD calculations upon which our 
work is based involve only two active flavours, the heavier third flavour is approximated 
by a renormalisation of the strong coupling constant. As a guide, recent work suggests 
that the best value of ex at may be 8 to 26 MeV larger than the value quoted above ||. 

One can notionally use QCD to directly calculate the value of cjv, but in practice 
the calculation has proven to be difficult. Early attempts to extract from the 
quark mass dependence of the nucleon mass (using Eq.(P) in quenched QCD with naive 
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extrapolations produced values in the range 15 to 25 MeV. Attention subsequently turned 
to determining by calculating the scalar matrix element of the nucleon (N\uu + dd\N) . 
There it was discovered that the sea-quark loops make a dominant contribution to 
[§,f|. These works, based on quenched QCD simulations found values in the 40 to 60 
MeV range, which are more compatible with the experimental values quoted earlier. 

On the other hand, the most recent estimate of a^, and the only one based on a 
two-flavour, dynamical-fermion lattice QCD calculation, comes from the SESAM collabo- 
ration. They obtain a value of 18 ± 5 MeV fl(| , through a direct calculation of the scalar 
matrix element (N\uu + dd\N) and the quark mass rh. 

The fact that neither (N\uu + dd\N), nor m is renormalisation group invariant in- 
troduces a major difficulty in calculating the sigma commutator in this approach. One 
must reconstruct the scale invariant result from the product of the scale dependent ma- 
trix element and the scale dependent quark masses. The latter are extremely difficult to 
determine precisely and are the chief sources of uncertainty. Furthermore, since lattice 
calculations are made at quite large pion masses, typically above 500 or 600 MeV, one 
needs to extrapolate, in the pion mass down to the physical value at 140 MeV. An impor- 
tant innovation adopted by Dong et al. was to extrapolate (N\uu + dd\N) using a form 
motivated by chiral symmetry, namely a + brh 2 . Regrettably, the value of b used was not 
constrained by chiral symmetry and higher order terms of the chiral expansion were not 
considered. Furthermore, since the work was based on a quenched calculation, the chiral 
behaviour implicit in the lattice results involves incorrect chiral physics JIT . 



3. THE CURRENT CALCULATION 

Our recent work [|12j was motivated by the improvements in computing power, together 



with the development of improved actions []13 1 , which have led to accurate calculations 



of the mass of the nucleon within full QCD (for two flavours) as a function of rh down 
to ~ 500 MeV. (Since m 2 is proportional to fh over the range studied we choose to 
display all results as a function of m. 2 .) We showed that provided that one has control 
over the extrapolation of this lattice data to the physical pion mass, one can calculate 
from <jn = m^dMpj / drn^ (which is equivalent to Eq. ([]]) where rh was used) at m v = 140 
MeV. This approach has the important advantage that one only needs to work with 
renormalization group invariant quantities. 

Chiral perturbation theory (%PT) predicts that the leading non-analytic (LNA) correc- 
tion to the self energy contribution to the nucleon mass is proportional to m\ (or m 3 / 2 ). 
It can be seen in Fig. that the preliminary point from CP-PACS |T3 at m 2 ~ 0.1 GeV 2 



does indeed suggest some curvature in this low mass region. These observations led the 
CP-PACS group to extrapolate their data with the simple, phenomenological form: 



M N = a + (3mi + 7m 3 , (2) 

rather than a naive linear form (7 = 0), as shown in Fig. [Tj. The corresponding fit to the 
combined data set, using Eq. (§), is shown as the short-dashed curve in Figs. [I] and We 
found that this fit gives o~n = 29.7 MeV. The difficulty with this purely phenomenological 



analysis was discussed in Ref. |L6[]. The problem is that a derivative is required when 
evaluating and the value of 7 found in the fit (-0.761 GeV -2 ) is almost an order of 




Figure 1. Nucleon mass versus m|. The solid 



data points are CP-PACS results [14 



the open points are UKQCD data |L5 



whilst 
. Both 

curves are fits using Eq. (FJ). The solid curve 
has 7 = 0, whilst the short-dashed curve has 
7 unconstrained. The vertical line indicates 
the physical pion mass. 



Figure 2. Data as labelled in Fig. [I]. The 
solid curve is a fit to Eq. (|3]) with a dipole 
form factor, the dashed curve is the same 
fit using a sharp cut-off form factor. The 
long-dash curve is a fit to Eq. fl3|) excluding 
the lowest data point. 



magnitude smaller than the model independent LNA coefficient, 7 LNA = —5.60 GeV~ 2 , 
indicated by ;\PT. 

Recently, an alternative approach was suggested in Ref. ||16|| . There it was realised 



that the pion loop diagrams shown in Fig. |] yield not only the most important non- 
analytic structure, but also give rise to the most significant variation in the nucleon mass 
as — > 0. This leads to the following extrapolation function for M^: 

M N = a + (3ml + cr NN {m^, A) + a NA (m n , A) , (3) 

where ct^n and ctjva are the self-energy contributions of Figs. 0(a) and 0(b), respectively, 
using a cut-off in momentum controlled by A. The full analytic expressions for u^n 



and o"ata are given in Ref. fig] . For our purposes it suffices that they have precisely the 
correct LNA and next-to-leading non-analytic behaviour required by chiral perturbation 
theory as — > 0. In addition, u^a contains the correct, square root branch point 
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Figure 3. One- loop pion induced self energy of the 
nucleon. 
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(M A -M N ) 2 }%) at the A-N mass difference, which is essential for extrapolations 
from above the A — Ntt threshold. 

Fitting Eq. (H) to the data, including the point near 0.1 GeV 2 , gives the dot-dash curve 
in Fig. 0. The corresponding value of is 54.6 MeV and the physical nucleon mass is 
870 MeV. Omitting the lowest data point from the fit yields the long-dash curve in Fig. || 
with o"tv = 65.8 MeV, demonstrating the need for lattice simulations of QCD at light 
quark masses. 



4. CONCLUSION 

The importance of the inclusion of the correct chiral behaviour is clearly seen by the fact 
that it increases the value of the sigma commutator from the 30 MeV of the unconstrained 
cubic fit to around 50 MeV. Nevertheless, it is a remarkable result that the present lattice 
data for two-flavour dynamical-fermion QCD, yields a stable [12[] and accurate answer for 
the sigma commutator, an answer which is already within the range of the experimental 
values. 

This work was supported by the Australian Research Council. 
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